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INTRODUCTION 


Integrated  optics  (10)  involves  the  technology  of  guiding  and 
manipulating  optical  waves  in  dielectric  waveguides,  the  generation  and 
detection  of  optical  waves,  and  the  coupling  of  optical  waves  into  and 
out  of  10  circuits.  The  simplest  form  of  optical  waveguides  is  a die- 
lectric slab.  By  studying  the  properties  of  slab  waveguides,  an  under- 
standing of  more  complicated  waveguides  is  possible. 

A planar  dielectric  slab  waveguide  is  shown  in  Figure  1.  The 
waveguide  region  has  an  index  of  refraction,  n, , with  surrounding  media 


having  indices  of  nQ  and  n2. 
must  be  met: 


ni  > n; 


£ nr 


For  waveguidlng,  the  following  criteria 

(1) 


FIGURE  1.  A Planar  Dielectric  Slab  Waveguide. 


1 Naval  Weapons  Center.  Integrated  Optics , by  L.  D.  Hutcheson. 
China  Lake,  Calif.,  NWC,  December  1975.  (NWC  TM  2592,  UNCLASSIFIED.) 
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If  n2  ■ n0,  then  it  is  called  a symmetric  waveguide.  If  n2  l*  nQ,  then 
there  exists  an  asymmetric  waveguide.  Of  course,  the  modal  solutions 
which  satisfy  the  boundary  conditions  are  simpler  for  the  symmetric  case 
than  for  the  asymmetric  case.  It  is  also  interesting  to  note  that  for 
the  symmetric  waveguide  there  is  no  cutoff  frequency2  which  implies  that 
the  lowest  order  mode  can  propagate  at  very  low  frequencies.  However, 
there  does  exist  a low  frequency  cutoff  for  all  the  modes  of  the  asym- 
metric waveguide. 

There  are  two  different  types  of  electromagnetic  modes  in  dielectric 
waveguides,  guided  modes  and  radiation  modes.  Figure  2 distinguishes 
between  these  two  different  types  of  modes.  To  obtain  total  Internal 
reflection,  the  incident  angle,  0j,  must  be  greater  than  the  critical 
angle,  0C,  for  each  interface.  For  interface  0,1 


For  interface  1,2 

(3) 


In  Figure  2(a)  the  light  passes  through  the  boundary  into  the  waveguide 
and  through  to  the  next  boundary;  this  is  called  a radiation  mode. 

Figure  2(b)  shows  another  radiation,  or  substrate,  mode  where  the  light 
is  incident  at  an  angle  larger  than  9C  for  the  upper  boundary  but  smaller 
than  "c  for  the  lower  boundary  with  n2  being  greater  than  nQ.  In  Figure 
2(c)  the  angle  0,  surpasses  both  critical  angles  of  each  Interface, 
which  results  in  total  internal  reflection.  Thus,  in  the  guided  mode 
the  light  is  trapped  in  the  guide. 

For  guided  modes  and  for  a certain  waveguide  thickness,  d,  and 
indices  no,  nj  and  n2,  light  will  propagate  in  the  guide  with  an  angle 
of  incidence  0 only  if,  after  two  successive  reflections,  the  wave 
front  is  again  in  phase  with  the  original  wave  front.  If  this  were  not 
the  case,  after  many  reflections,  somewhere  down  the  waveguide  wave 
fronts  with  a range  of  phases  between  0 and  2”  would  add  to  zero  ampli- 
tude, which  implies  the  wave  would  not  propagate.  This  results  in  the 
requirement  that  at  an  arbitrary  point  the  phase  of  one  wave  front 
obtained  from  another  by  two  successive  reflections  must  equal  the  phase 
of  the  second  wave  front  at  the  same  point,  or  must  differ  by  a multiple 


2 D.  Marcuse.  Theory  of  Dielectric  Optical  Waveguides.  New  York, 
Academic  Press,  1974. 

. 3 R.  A.  Andrews.  Optical  Waveguides  and  Integrated  Optics  Tech- 
nology, Naval  Research  Laboratory  Report  7291  (August  1971),  UNCLASSIFIED. 
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of  2n.  This  requirement  permits  propagation  for  only  discrete  values  of 
li.  This  means  that  only  a finite  number  of  discrete  guided  modes  exists. 
This  is  quite  different  from  the  radiation  modes  in  that  there  exists  an 
infinite  continuum  of  unguided  radiation  modes. 

This  report  will  present  the  waveguiding  properties  of  optical 
dielectric  waveguides  using  both  the  ray  optics  and  wave  optics  tech- 
niques. Then  these  results  will  be  applied  to  investigate  the  losses 
in  dielectric  waveguides.  In  particular,  losses  due  to  short  bending 
radii  of  the  dielectric  waveguides  will  be  presented. 


THEORY  OP  DIELECTRIC  WAVEGUIDES 


MAXWELL'S  EQUATIONS 


Maxwell's  equations  are  explained  in  detail  in  many  textbooks4'5 
on  electromagnetism,  so  here  they  will  only  be  stated  and  used  ^s  needed. 
Maxwell's  equations  relate  the  magnetic  field  intensity  vector  H and  the 
magnetic  flux  density  vector  B to  the^electric  field  intensity  vector  E 
and  the  e'  metric  displacement  vector  D: 


-> 

V 


> 

H 


3D 

3t 


(4) 


and 


-► 

7 x E 


3B 

3t 


(5) 


where  3/3t  is  the  partial  derivative  with  respect  to  time  and  V has  the 
vectorial  components  of  [3/3x,  3/3y,  3/3z).  (Maxwell's  equations  will 
be  used  in  the  context  of  light  wave  interactions.)  Consequently,  the 
current  terra  has  been  neglected  in  the  above  equation.  Also,  it  will 
be  assumed  that  there  are  no  electric  charges  which  means  that  D satis- 
fies the  relation 


V • 6 - 0 


(6) 


Similarly,  B satisfies  the  relation 

> > 

V • B - 0 


(7) 


In  addition  to  Maxwell's  equations  there  are  equations  which  relate 
the  two  electric  vectors  and  the  two  magnetic  vectors  which  are,  in  gen- 
eral, material  dependent: 


4 C.  H.  Dumey  and  C.  C.  Johnson.  Introduction  to  Modem  Electro- 
magnetics} . New  York,  McGraw-Hill,  1969. 

5 J.  A.  Stratton.  Electromagnetic  Theory.  New  York,  McGraw-Hill, 

1941. 


4 


NWC  TM  3136 

B - e • E (8) 

and 

B = v • H , (9) 

where  £ and  ij  are  tensors  representing  the  material  permittivity  and 
permeability,  respectively.  If  £ and  il  are  functions  of  positions  in 
the  medium  then  the  medium  is  called  inhomogeneous.  If  there  is  no 
position  dependence,  the  medium  is  called  homogeneous.  On  the  other 
hand,  if  £ and  M are  tensors,  the  medium  is  termed  anisotropic  while  an 
isotropic  medium  is  one  in  which  £ and  u are  scalar  quantities.  There- 
fore, the  most  general  type  of  medium  is  anisotropic  and  inhomogeneous. 
However,  for  many  cases  of  practical  interest,  the  medium  can  be  approx- 
imated to  be  isotropic  and  homogeneous,  which  is  what  will  be  assumed  in 
the  following  analysis. 

To  obtain  the  useful  wave  equation  one  must  manipulate  Maxwell's 
relations  until  an  equation  results  which  depends  only  on  ft  or  ft.  To 
make  the  analysis  simpler,  without  losing  generality,  one  can  assume 
one  of  the  guided  dimensions  to  be  much  larger  than  the  other  guided 
dimension.  In  the  waveguide  shown  in  Figure  1,  the  y-dimension  is  made 
much  larger  than  the  x-direction  which  yields  a planar  waveguide  and  in 
Maxwell's  equations  one  can  take 


The  planar  waveguide  is  the  same  solution  as  the  channel  waveguide  with 
one  dimension  much  larger  than  the  other. This  actually  does  not 
restrict  the  generality  of  the  modal  analysis  since  the  coordinate  sys- 
tem can  always  be  rotated  until  this  condition  is  satisfied  for  any 
given  mode. 

Before  one  can  fully  understand  the  interactions  of  guided  waves  a 
knowledge  of  the  properties  of  guided  modes  must  be  obtained.  A concise 
definition  of  the  term  "mode"  is  difficult  to  make.  A mode  can  be 
regarded  as  an  eigensolution  of  Maxwell's  equations  belonging  to  a 
particular  eigenvalue  and  satisfying  all  the  boundary  conditions  of  the 
problem. ' The  modes  of  a dielectric  slab  waveguide  are  classified  as 
transverse  electric  (TE)  and  transverse  magnetic  (TM)  modes.  A TE  mode 


( E.  A.  J.  Marcatlli.  "Dielectric  Rectangular  Waveguide  and 
Directional  Coupler  for  Integrated  Optics,"  BELL  SYST  TECH  J,  Vol.  48 
(September  1969),  pp.  2071-2102. 

1 D.  Marcuse.  Light  Transmission  Optics.  New  York,  Van  Nostrand 
Reinhold,  1972. 
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means  that  the  electric  field  component  in  the  direction  of  wave  propa- 
gation does  not  exist.  A TM  mode  has  no  component  of  the  magnetic  field 
in  the  longitudinal  direction.  Since  the  TE  and  TM  modes  have  different 
components  they  must  be  studied  separately.  • 

TE-GUIDED  MODES  * 


Transverse  electric  modes  only  have  three  field  components  which 
are  Ey,  Hx,  and  Hz.  In  a linear  and  isotropic  material  Maxwell's  Equa- 
tions 4 and  5 become,  after  using  relations  8 and  9, 


->  -►  o 3E 

V x H = e0n  — 

(id 

and 

V x E = -pfl  , 

(12) 

with 

7 = X 3^  + y 37  + 2 3l 

* 

(13) 

where  x,  y and  z are  unit  vectors  in  the  x,  y and  z directions,  respec- 
tively, and  n is  the  index  of  refraction  of  the  medium.  Assuming  a 
field  of  the  form 

i(wt  - (Jz) 

e » 

(14) 

where  8 is  the  propagation  constant  along  the  axis  of  waveguide  and 
using  relation  10,  Maxwell’s  equations  yield 

3H 

-l6“,  - “ST  ' ‘"‘o"2*. 

9 

(15) 

i6Ey  = -iu>p0Hx  , 

(16) 

3E 

if  • -l,*oHs  • 

(17) 

Putting  H and  H in  terms  of  E t 
x z y 

we  get 
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Substituting  Equations  18  and  19  Into  15  yields  the  familiar  wave  equa- 
tion In  terms  of  Ey, 

;>?e 

£+  (k2n2  - 82)  E = 0 , (20) 

3x2  y 

where  k2  = W2u  qC.  Therefore,  by  solving  the  one-dimensional-wave  Equa- 
tion 20  for  Ey,  the  magnetic  field  components,  and  Hz,  are  obtained 
directly  from  Equations  18  and  19. 

Before  trying  to  solve  the  wave  equation  in  each  of  the  three 
regions,  a physical  picture  of  the  modes  can  be  obtained  by  seeing  what 
happens  when  the  propagation  constant,  8,  at  some  fixed  frequencv,  takes 
on  different  values.  Refer  to  Figure  1 to  see  the  relationship  between 
the  three  regions  and  the  indices.  Also,  the  waveguide  is  assumed  to  be 
asymmetric  with  n(  > n2  > nQ.  Therefore,  the  wave  equation  is  different 
for  each  of  the  three  regions  since  n in  Equation  20  takes  on  nQ,  nj 
and  n2*  Clearly,  if  6 > knj , which  is  the  largest  0 possible,  and8 
since  1/E  32E/3x2  > 0,  then  E must  be  exponential  in  all  three  regions. 
The  boundary  conditions  at  both  interfaces  must  be  matched  which  means 
that  the  field  increases  to  infinity  exponentially  away  from  the  wave- 
guide. This  says  that  for  this  8 the  solution  is  not  physically  realiz- 
able and  is  not  a real  wave. 

To  obtain  a guided  mode  solution  the ^propagation  constant  must 
satisfy  knQ,  kn2  < 8 < knj  which  yields  a sinusoidal  solution  in  region 
one  and  an  exponential  solution  in  the  other  two  regions.  This,  of 
course,  can  be  achieved  only  if  nj  > n2,  nQt  which  further  supports  the 
criterion  that  the  waveguiding  layer  must  have  the  largest  index  as  was 
shown  in  Equation  1.  The  radiation  modes  exist  for  knQ  < 8 < kn2,  which 
is  the  substrate  radiation  mode.  Also,  for  0 < 6 < knQ  radiation  modes 
exist  in  both  surrounding  media  since  the  field  is  sinusoidal  in  all 
three  regions. 

When  solving  the  wave  equation,  the  solutions  must  satisfy  the 
boundary  conditions  at  the  two  dielectric  interfaces  at  x ■ 0 and  x “ 

-d.  The  boundary  conditions  are  such  that  the  tangential  E and  H fields 
must  be  continuous  at  the  dielectric  interfaces.  This  means  that  both 
Ey  and  Hz  must  be  continuous  at  x = 0 and  x * -d.  In  addition  to  this 
requirement,  to  obtain  a physically  realizable  solution  the  electric 
field,  Ey,  must  be  zero  at  x = ±°°.  Since  the  term  k2n2  - 82  takes  on 
different  values  in  each  of  the  three  regions,  it  will  be  simpler  to 
define  three  new  variables: 


A.  Yariv.  Quantum  Electronics,  2nd  Edition.  New  York,  John 
Wiley  and  Sons,  Inc.,  1975. 
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«2  = k2n2  - B2 
Y2  = 82  - k2n| 


62  = g2  - k2n2 


These  definitions  provide  for  <*,  Y and  6 being  real  and  the  solution  to 
the  wave  equation  can  be  assumed  to  be  of  the  form 


Ae-dx 

0 < x 

Beiax  + Ce~iotx 

-d  < x £ 0 

DeY(x+d) 

-d  i X 

The  constants  A,  B,  C and  D are  determined  through  the  use  of  the 
boundary  conditions  at  x = 0 and  x = -d,  upon  which  applying  the  conti- 
nuity of  Ey  across  the  boundary  yields 

A = B + C (25) 


D - B.'1Ud  + CelQd  . (2 t 

and  since  Hz  must  also  have  its  tangential  components  equal  across  the 
boundaries  the  following  additional  equations  are  obtained: 

-6A  = iuB  - iaC  (27 


6D  - iuBe~lud  - iaCeiad 


Now  we  have  four  equations  and  the  four  unknown  constants  (A,  B,  C and 
D).  Solving  the  four  equations  by  the  method  of  determinants  results 
in  the  following  relation: 

TYT  . r&n 


ud  ■>  arc  tan  Jj^-]  + arc  tan  J 


The  fact  that  Equation  29  is  periodic  with  a period  it,  since  the  tan  is 
periodic  in  n,  says  that  relation  29  should  be  written 

ud  ■ arc  tan["^  + arc  tanCa3  + ■*  » (30) 

where  m - 0, 1,2,3,  ...  is  the  mode  number. 
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Substituting  the  expressions  for  Ey  (Equation  24)  into  the  wave 


equation  (Equation  20),  yields 

<52  + k;’n2  - 62  = 0 , (31) 

-«2  + k2n2  - 62  = 0 , (32) 

. Y2  + k2n2  - B2  = 0 , (33) 

where  k **  2n/X0»  and  eliminating  6 from  the  equations  we  get, 

62  + a2  = k2 (n2  - n2)  , (34) 

Y2  + a2  = k2 (n2  - n2)  . (35) 


These  two  equations,  together  with  Equation  30,  are  the  complete  set  of 
equations  which  determine  the  TE  modes.  The  above  equations  can  be 
applied  to  either  an  asymmetric  or  a symmetric  guide.  First,  the  equa- 
tions will  be  solved  for  the  asymmetric  case.  For  a typical  asymmetric 
guide,  the  refractive  indices  nQ,  nj  and  n2  satisfy  the  relationship 
that  nj  >>  nQ  and  n1  s;  n2.  This  would  be  true  in  most  cases  when  the 
top  boundary  layer  is  air  (nQ  = 1).  The  simplest  way  to  solve  these 
equations  is  to  solve  them  graphically.  Relations  34  and  35  are  equa- 
tions of  a circle  with  radii  of 

= (62  + a2)*5  (36) 

and 

r2  = (y2  + a2)**  , (37) 

and  they  are  shown  in  Figure  3.  For  the  indices  of  nQ.  = 1,  nj  * 1.55 
and  n2  = 1.52,  the  circle  describing  62  + a2  would  be  approximately  15 
times  larger  than  y2  + ot2.  As  can  be  seen  from  Figure  3,  the  range  of 
values  which  a and  y can  take  on  are  the  same  and  are  much  larger  than 
range  of  6.  However,  the  magnitude  of  6 is  much  larger  than  the  magni- 
tudes of  u and  y.  This  will  simplify  Equation  30  by  using  6/a  » 1 
which  makes  arc  tan  6/a  » it/2  and  leads  to 

ad  **  arc  tan{^aJ  + + inn  , (38) 

which  can  be  rewritten  in  the  form 

yd  * ad  tan  (ad  - ir/2  - on)  ■ -ad  cot  ad  . . (39) 


9 


NWC  TM  3136 


FIGURE  3-  Diagram  Which  Illustrates  the  Relative 
Magnitude  of  a,  6 and  y. 

Equation  39  and 

(Yd)2  + (ud)2  = (kd)2(n2  - n2)  (40) 

are  each  plotted  in  Figure  4 and  the  intersection  of  the  two  curves  are 
the  solutions.  As  discussed  earlier,  y must  be  positive  to  have  guided 
modes;  consequently,  only  the  first  quadrant  is  plotted  for  the  circles. 
The  intersections  where  y is  negative  correspond  to  the  even  modes; 
therefore,  only  odd  modes  exist  for  the  asymmetric  waveguide.  It  can 
be  seen  that  asymmetric  dielectric  waveguides  have  a low  frequency 
(long  wavelength)  cutoff  below  which  energy  does  not  propagate.  For  a 
particular  propagating  mode  the  values  of  y and  a are  determined  by  the 
circle  of  radius  R where 

R - kd(n2  - n2)*5  , (41) 

and  for  the  mth  propagating  mode  the  radius  of  the  circle  must  be  between 
(m  + 1)"  i kd(n‘j  - n2)  i (2m  + l)n/2  . (42) 
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FIGURE  4.  Graphical  Solutions  of  the  Transcendental 
Equations  for  the  Asymmetric  Waveguide. 


The  condition  of  index  change  An  = 
becomes 


An  ** 


(2m  ~ l)2  f 02 

32n2  d 


9 


n2  for  cutoff  of  the  m1*1  mode 

(43) 


where  \ is  the  free  space  wavelength  and  d is  the  waveguide  thickness. 
The  number  of  modes  that  propagate  are  finite  and  are  dependent  on  the 
parameters  'nj , n2,  X and  d which  are  controllable.  For  example,  keep- 
ing X/d  constant,  an  index  change  by  a factor  of  nine  is  required  for 
the  second  mode  to  propagate  over  that  required  for  the  lowest  order 
mode. 


If  the  dielectric  waveguide  is  symmetric,  i.e.,  n2  ■ nQ,  then  some 
surprisingly  different  results  are  obtained.  The  analysis  is  simplified 
if  the  modes  arc  separated  into  even  and  odd  modes  due  to  the  symmetry. 
Considering  even  modes  first,  from  Equation  24, Ey  in  the  waveguide 
becomes  (subscript  e is  for  even) 


* 
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E = B cos  ux 

y e 

and  from  Equation  19, 
iu 


H_ 


B sin  ax 
wu  0 e 


Ey  In  the  bottom  layer  is 


E = D cos  ad  e 

y e 


Y (x+d) 


and 


iY 


z uiu0  e 


D cos  ad  eY(x+d) 


(44) 


(45) 


(46) 


(47) 


Applying  the  boundary  condition  at  x = -d,  that  E must  be  continuous 
results  in  ^ 


and  applying  continuous  at  x * -d  yields 
Yd  = ad  tan  ad 


(48) 


(49) 


The  amplitude  Be  can  be  determined  by ^calculating  the  power  flowing 
in  the  mode  along  the  waveguide  (z-direction)  using  Poynting’s  theorem 

P - J f~(E  x H*)'z  dx  . (50) 

«oo 

Since  the  above  involves  a cross  product  and  the  only  component  of  a TE 
mode  is  Ey,  then  the  magnetic  field  component  needed  is  Hx;  therefore. 


? m - — f E H*  dx 

2 _ y * 


and  from  Equation  18, 

6 


H 


which  gives 


E 


y 


(51) 


(52) 


P - 


wu . 


r iej2  ^ . 


(53) 


Substituting- Equations  44  and  46  into  Equation  53,  and  using  Equation 
49,  yields' 
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r 2w%p  i5* 

[fid  + 6/yJ 


Similarly,  the  odd  modes  are  solved  in  the  same  manner  using 
E = C„  sin  ax 

y 0 

inside  the  waveguide  with 

H = r~~  Cn  cos  ax 

Z wp  o 0 

In  the  bottom  layer. 


E = -En  sin  ad  e 
y 0 


Y(x+d) 


,,  -iY  4 , Y (x+d) 

H = - — E0  sin  ad  e , 


where  the  minus  sign  comes  from  the  fact  that  it  is  an  odd  mode.  Using 
the  continuity  of  the  fields  across  the  boundary  yields 

C0  = E0  (59) 


Yd  « -ad  cot  ad 

The  power  flow  calculation  yields 

r ,p  "i** 


2u»u0P  p 

6d  + 6/yJ 


One  has  to  be  careful  when  using  Equations  54  and  61  since  they  look 
the  same.  When  using  the  relation  for  Bg,  the  constants  6 and  Y are 
evaluated  from  relation  49;  whereas,  relation  60  is  used  to  evaluate  6 
and  Y for  CQ. 

Figure  5 shows  a plot  of  relations  49  and  60,  and  relation  40  is 
used  in  the  same  manner  to  find  the  graphical  solutions  as  was  done  for 
the  asymmetric  waveguide.  -As  can  be  seen  from  the  solutions  of  th* 
transcendental  equations  for  the  symmetric  waveguide,  there  is  no  low 
frequency  cutoff  for  the  lowest  order  propagating  mode.  The  frequency 
can  be  made  arbitrarily  low  (smaller  radius  circle)  and  a solution 
still  exists  (see  Figure  5). 
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(EVEN  MODES)  (ODD  MODES) 


FIGURE  5.  Graphical  Solutions  of  the  Transcendental 
Equations  for  the  Symmetric  Waveguide. 


Another  Important  aspect  of  guided  wave  optics  is  the  dispersion 
in  the  waveguides,  i.e.,  how  the  propagation  constant,  6,  varies  with 
frequency.  The  dispersion  relationship  is  obtained  by  substituting 
the  equations  for  a,  Y and  6 (Equations  21-23)  into  Equation  30,  which 
yields 


l(k2nj  - 62)d2l*4  - arctan 


(B2  - k^2)*5 


+ arctan 


<k2n2  - e2)1* 

V -Vo2)* 

(k2n2  - 


(62) 


+ mu 


for  TE  inodes  of  an  asymmetric  waveguide.  Also,  remember  for  guided 
modes,  the  propagation  constant  obeys 


knfl  < kn2  < 8 < knt 


(63) 
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where  k = <*>/c.  Therefore,  as  can  be  seen  In  Figure  6,  6 is  bounded  by 
w/c  n(  above  which  no  solutions  exist,  and  w/c  n2  below  which  exist 
radiation  modes.  Figure  6 also  shows  what  three  typical  dispersion 
equations  look  like,  and  Figure  7 shows  computer  solutions  obtained  for 
an  asymmetric  waveguide  with  nQ  ■ 1,  n,  « 2.23,  n2  = 1.393  and  d «•  3.34 
pm.  **  Physically,  the  shape  of  the  dispersion  curves  is  explained  as 
follows.  Increasing  the  frequency,  the  depth  of  penetration  of  the 
exponential  tails  of  the  fields  into  the  surrounding  media  decreases, 
which  means  that  the  electric  fields  are  more  confined  to  the  guiding 
layer.  In  the  limit  as  w -*■  “ the  fields  are  totally  confined  to  the 
waveguide  layer.  The  slope  of  the  upper  boundary  is  the  velocity  of 
light  in  the  substrate  with  index  n2.  When  w/c  approaches  kn2  the 
constant  Y approaches  zero  and  this  causes  the  exponential  tail  of  the 
field  to  extend  infinitely  into  the  substrate. 


FIGURE  6.  Typical  w Versus  B Diagram  Showing  the  Dispersion 
Relationship  of  Three  Discrete  Guided  Modes  for  Dielectric 
Slab  Waveguides. 


0 M.  K.  Barnoski.  Introduction  to  Integrated  Optics.  New  York, 
Plenum  Press,  1974. 
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FIGURE  7.  Dispersion  Curves  of  Ten  Modes  of  an 
Asymmetric  Waveguide  With  d = 3.34  pm,  nQ  * 1, 
nt  ■ 2.23  and  n2  ■ 1.393.  Natural  dispersion  is 
neglected. 

Beyond  the  limits  of  6 * w/c  n2  and  B = w/c  nj , there  is  a continuous 
spectrum  of  radiation  modes.  In  this  region  6 is  imaginary  and  these 
modes  are  called  evanescent.  Figure  8 shows  typical  electric  field 
patterns  for.  the  three  lowest  order  modes  of  an  asymmetric  waveguide. 


This  concludes  the  discussion  for  the'TE-guided  modes  of  a dielec- 
tric waveguide.  The  analysis  for  TM  modes  is  very  similar  except  t*>e 
field  components  are  Ex,  Ez  and  Hy  The  TM  modes  are  obtained  by  set- 
ting Hz  * 0.  The-  resulting  transcendental  equations  are  much  the  same. 
The  equations  of  the  circles  are  identical  to  those  for  the  TE  modes 
(i.e.t  relations  34  and  35).  The  difference  between  TE  and  TM  modes  is 
in  Equation  30,  where  the  ratio  y/a  in  the  TE  case  becomes  (ni/n2)2  y/o 
for  TM  modes  and  the  ratio  6/a  in  the  TE  case  becomes  (n^n^2  6/a  for 
TM  modes. 
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X 0 -d 


FIGURE  8.  Typical  Electric  Field  Intensity 
Patterns  for  Three  Lowest  Order  Modes. 


‘ Now  we  have  a useful  understanding  of  guided  wave  optics.  How  the 
parameters  of  the  waveguide  can  be  changed  to  affect  the  propagating 
modes  is  understood.  This  treatment  can  be  extended  to  two  dimensions 
(i.e.,  rectangular  dielectric  waveguides)  but  the  results  are  much  the 
same  except  guiding  takes  place  in  both  of  the  transverse  directions. 

In  deriving  the  wave  equation,  3/3 y can  no  longer  be  taken  as  zero.  At 
any  rate,  now  we  can  consider  other  important  topics  by  using  our  knowl- 
edge of  the  waveguide  modes.  The  next  section  will  consider  losses  in 
dielectric  waveguides  due  to  bending  of  the  waveguides.  These  losses 
are  important  in  establishing  the  maximum  packing  density  of  integrated 
optical  circuits  (IOCs). 


BENDING  LOSSES 

There  are  many  factors  which  limit  the  performance  of  waveguides 
in  IOCs.  Two  of  the  more  important  ones  are  scattering  losses  which 
are  due  to  waveguide  edge  roughness,  and  bending  losses  which  are  due 
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to  a short  radius  of  curvature  of  the  waveguide.  Bending  losses  is  not 
necessarily  the  more  important  of  the  two  losses  in  describing  the 
performance  limitations  of  the  guides.  However,  if  one  is  to  be  con- 
cerned with  the  lower  limit  of  the  number  of  functions  per  unit  area  to 
be  performed  on  an  IOC,  then  the  primary  topic  to  be  considered  is  bend- 
ing losses. 

No  matter  how  large  or  small  the  radius  of  curvature  of  the  bend, 
the  dielectric  waveguides  will  always  lose  power  by  radiation  around 
the  bend.  The  loss  may  be  quite  small  if  the  radius  of  curvature  is 
sufficiently  large.  However,  losses  can  become  quite  large  if  the 
radius  of  curvature  is  small.  The  bend  losses  come  from  two  effects: 
First,  the  shape  of  the  mode  varies  with  the  radius  of  curvature,  and 
second,  the  energy  is  radiated  to  the  surrounding  media  as  it  travels 
around  the  bend.  The  bending  loss  is  a very  fast  function  of  the  radius 
of  curvature.  In  other  words,  the  loss  due  to  some  value  of  radius  of 
curvature  can  be  quite  small,  but  a very  small  change  in  the  bend  radius 
will  cause  a sharp  increase  in  the  loss. 

First,  let  us  get  a physical  picture  of  the  causes  of  the  bending 
losses.  Figure  9 shows  a waveguide  with  a radius  of  curvature  R.  A 
portion  of  the  energy  of  the  dielectric  waveguides  travels  outside  the 
waveguide  in  the  form  of  an  evanescent  wave.  The  evanescent  field  out- 
side the  waveguide  is  decaying  exponentially  in  the  transverse  direc- 
tion to  the  waveguide  axis.  As  the  waveguide  is  bent  the  planes  of 
constant  phase  are  as  shown  in  Figure  9.  The  planes  of  constant  phase 
are  separated  further  apart  the  farther  one  goes  from  the  center  of 
curvature.  Therefore,  the  farther  from  the  center  of  curvature  the 
wave  is,  the  faster  it  must  travel  so  that  it  can  still  be  in  phase 
with  the  guided  wave.  Each  medium  has  its  own  characteristic  phase 
velocity  which  is  determined  by  the  index.  The  surrounding  media  has 
a larger  phase  velocity  than  does  the  waveguide  since  the  index  is 
smaller.  At  some  critical  distance  from  the  waveguide,  the  evanescent 
wave  must  travel  faster  than  its  characteristic  phase  velocity  to  remain 
in  phase  with  the  rest  of  the  wave.  Since  this  is  not  a physically 
realizable  situation  the  wave  detaches  itself  from  the  rest  of  the  wave 
and  radiates  into  space.  The  sharper  the  bend  of  the  waveguide,  the 
closer  the  critical  distance  is  to  the  waveguide.  Since  the  evanescent 
wave  is  an  exponentially  decaying  wave  the  field  increases  exponentially 
towards  the  waveguide  and,  as  will  be  shown,  the  critical  distance  varies 
linearly  with  the  radius  of  curvature.  Therefore,  only  a small  change 
In  the  radius  can  make  an  exponential  change  in  the  field  that  is  radi- 
ated away  from  the  waveguide. 

Figure  10  shows  a propagating  mode  along  with  its  associated 
evanescent  wave  in  an  expanded  view  of  a curved  waveguide.  The  angle  9 
Is  the  angle  between  adjacent  planes  of  constant  phase.  At  the  critical 
distance  (R  + r),  the  distance  that  the  evanescent  wave  must  travel  from 
one  plane  of  constant  phase  to  the  next  is  (R  + r)6.  The  distance  that 
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the  guided  wave  must  travel  Is  RQ.  However,  they  must  travel  these 
different  distances  in  the  same  amount  of  time.  The  velocity  inside 
the  waveguide  is  w/knj  and  the  maximum  velocity  permissible  in  the  outer 
media  is  u>/kn?.  By  setting  this  maximum  velocity  at  the  critical  dis- 
tance and  equating  the  times  for  the  waves  to  travel  the  different  dis- 
tances, we  get 

(R  + r)  0kn2  RekOj 


and  solving  for  r,  the  critical  distance,  we  get 


The  attenuation  loss  for  slab  waveguides  can  be  calculated  using 
the  velocity  approach  described  above  and  first  developed  by  Marcatili 
and  Miller.  0 From  the  discussion  of  slab  waveguides  in  the  previous 
section,  the  transverse  field  distribution  can  be  described  as 

E * cos(kxx)  for  ~ £ x < j 

and 

E *=  cos(kxd)e_[lxl  " d]/n  for  ‘ |x|  i f , (66) 

• i 

where  the  coordinate  system  is  now  in  the  middle  of  the  waveguide,  n~l 
is  the  transverse  decay  coefficient  outside  the  guide,  and  the  waveguide 
is  symmetric.  Obviously,  to  discuss  losses  due  to  bends,  as  described 
by  Figure  10,  the  waveguide  must  have  guidance  in  two  dimensions.  How- 
ever, in  this  analysis  the  losses  in  the  plane  of  the  waveguide  will  be 
considered.  Later  in  this  section  there  will  be  a discussion  of  two- 
dimensional  losses. 


Now  split  the  energy  carried  by  the  field  into  two  parts:  (1)  that 
energy  beyond  the  critical  distance  r(£r)  and  (2)  the  total  energy  in 
the  waveguide  (Sw)»  This  yields  the  two  energies10 

£r  - n/2  cos2(kxd/2]e“2fr  ' b/2l/n  (67) 

and 

£w  - n/2  + l/2kx  sin  kxd  + n cos2[kxd/2]  . (68) 


10  E.  A.‘  J.  Marcatili  and  S.  E.  Miller.  "Improved  Relations 
Describing  Directional  Control  in  Electromagnetic  Waveguidance,"  BELL 
SYST  TECH  J,  Vol.  48  (September  1969),  pp.  2161-88. 
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Now  let  PQ  be  the  propagating  power  and  a the  attenuation  constant 
where 


1 dP 


2PC  dz 


(69) 


* “ tx 


and  the  use  of  the  relations  for  the  energy  (Equations  67  and  68)  yields 

(71) 


[d  + 2ncos(dxd/2)]2 


2X 


The  attenuation  constant  a In  terms  of  the  energy  is  given  by 


1 fir 


U £T 


(72) 


Substituting  in  the  relations  for  r,  l,  and  yields 
-BR 


a ■ A e 


(73) 


where 


[Xn/2]cos2 (k  d)eT1^5^ 


[n/2  + (l/2k  )sink  d + ncos2(k  d/2)][d  + 2ncos(k  d/2)] 


(74) 


and 


2fni  - °n 

{ n2  J 


(75) 


These  results  further  show  the  commanding  influence  that  the  bend- 
ing radius  has  on  the  loss.  Table  1 shows  some  results*®  using  the 
above  analysis  with  the  waveguide  index  of  refraction  ■ 1.5.  Figure  11 
shows  a plot  of  the  loss  as  a function  of  R for  case  2.  Note  that  a 
40Z  decrease  in  R brings  about  two  orders  of  magnitude  increase  in  the 
l0S8  . 


Marcatill  has  made  a complete  analysis  of  the  losses  in  rectangular 
dielectric  waveguides.**  In  his  analysis  he  used  a modal  analysis  as 


11  E.  A.  J.  Marcatill,  "Bends  in  Optical  Dielectric  Guides,"  BELL 
SYST  TECH  J,  Vol.  48  (September  1969),  pp.  2103-32. 
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Let  a be  the  transverse  field  width,  then  t,  the  distance  that  the  energy 
remains  collimated  in  an  infinite  medium,  is  given  by 
.2 

» (70) 
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opposed  to  the  velocity  approach  described  above.  Since  his  paper  has 
such  a complete  analysis,  and  his  results  are  very  similar  to  those  just 
obtained,  this  discussion  will  just  present  some  of  the  important  results 
not  covered  in  the  velocity  analysis. 

In  the  modal  analysis,  the  effect  of  guide  curvature  is  to  change 
the  real  part  of  the  propagation  constant,  introduce  an  imaginary  part 
of  the  propagation  constant  corresponding  to  the  loss,  and  alter  the 
field  distribution.  This  also  leads  to  the  radiation  loss  being  related 
to  the  bend  radius  exponentially.  Figure  12  shows  the  field  distribution 
as  a function  of  guide  width,  d.  It  also  shows  that  the  field  distribu- 
tion perpendicular  to  the  plane  of  curvature  is.  symmetric  and  unchanged 
due  to  bending  of  the  waveguide. 


FIGURE  12.  Field  Distribution  as  a Function 
of  Guide  Width  d.  All  three  waveguides  have 
the  same  radius  of  curvature. 
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DISCUSSION 


It  has  been  shown  that  losses  in  dielectric  optical  waveguides  due 
to  bends  in  the  waveguides  can  be  quite  large.  When  considering  packag- 
ing density  of  integrated  optical  circuits,  bending  losses  definitely 
need  to  be  taken  into  account.  The  losses  are  an  exponential  function 
of  the  bending  radius. 

Further  investigation  of  the  bending  losses  needs  to  be  performed. 
In  particular,  experimental  studies  are  needed.  A number  of  theoretical 
aspects  of  bending  losses  have  been  reported;  however,  no  experimental 
investigations  have  yet  been  reported.  The  next  step  of  this  program  is 
to  investigate  the  experimental  aspects  of  bending  losses.  A program  is 
now  underway  to  make  diffused  waveguides  in  LiNbOj  and  to  experimentally 
characterize  the  waveguides  and  compare  the  results  with  theory.  Once 
this  has  been  accomplished,  curved  waveguides  will  be  made  and  charac- 
terized, with  an  emphasis  placed  on  the  dependence  of  bending  loss  on 
bending  radius. 


